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, We have investigated the phase behaviour of a symmetrical binary fluid mixture for the situation 

^1^' where the chemical potentials /ii and /i2 of the two species differ. Attention is focused on the set of 

interparticle interaction strengths for which, when /xi = /i2, the phase diagram exhibits both a liquid- 
, vapor critical point and a tricritical point. The corresponding phase behaviour for the case fii 7^ fi2 

. is investigated via integral-equation theory calculations within the mean spherical approximation 

(MSA), and grand canonical Monte Carlo (GCMC) simulations. We find that two possible subtypes 
of phase behaviour can occur, these being distinguished by the relationship between the critical lines 
in the full phase diagram in the space of temperature, density, and concentration. We present the 
i-C ' detailed form of the phase diagram for both subtypes and compare with the results from GCMC 

simulations, finding good overall agreement. The scenario via which one subtype evolves into the 
other, is also studied, revealing interesting features. 
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I. INTRODUCTION 



I As one proceeds from a simple single component fluid to a binary mixture, the complexity of the phase behaviour 
• increases dramatically. The range of possible behavior was first summarized on a qualitative level more than 25 
I [ years ago in the seminal study of van Konynenburg and Scott Using a simple mean- field (MF) ansatz for the 
''O • equation-of-state, these authors found that in binary mixtures the additional thermodynamic degree of freedom in 
Ph ■ the form of the concentrations of the two species, ci and C2, with c = ci — 1 ~ C2 (or, alternatively, the difference in 
the chemical potentials of the two species, = /ii — /i2), considerably widens the spectrum of possible features 
Among these are upper and lower critical points, tricritical points, critical end points (CEPs), quadruple (i.e., four 
phase) points, lines of critical points, or triple lines. The details of which features actually emerge in practice for a 
given model, and in which order (i.e., the phase diagram topology), were found to depend in a non-trivial manner 
on the precise combination of interparticle interaction parameters, i.e., three energy parameters and three length 
scales aij (with i,j — 1, 2), together with parameters that characterize the range of the attractive interactions. 
' Today, more than a quarter of a century later, both liquid state theories _3J and computer simulations 4, 5] offer 
Q\ . standards of accuracy that greatly surpass those obtainable from MF-based theories. It therefore seems appropriate 
' to reconsider the schematic overview of phase diagram topologies provided by van Konynenburg and Scott p|, at 
, a more quantitative level. The present contribution aims to take a step in this direction. However, owing to the 
high-dimensional space of interaction parameters, one is necessarily restricted to a particular model. In the interests 
of numerical tractability, we have chosen to consider a symmetrical binary mixture. In such a system, the interactions 
between particles of like species are equal, i.e., $ii(r) — ^22{r) — ^{r), while the interactions between unlike species 
are related via $12 (r) = S^{r). Thus for a prescribed form of interaction potential, <I>(r), the phase behaviour is 
controlled by the single parameter S. Since for the case S > 1, such a system exhibits only a simple liquid- vapour 
(LV) transition, we focus here on the more interesting case of 5 < 1, where the competition between a LV on the one 
Q hand, and a demixing transition (into phases of generally different concentrations), gives rise to a rather rich phase 
O ' behaviour. For the reasons explained below, we will focus on a mixture where $(r) takes the form of a hard-core plus 
I attractive Yukawa (HCY) interaction; nevertheless we anticipate that any other similar forms of potential (such as 
" '~j ■ Lennard-Jones or square-well potentials) will lead to results qualitatively similar to those we present. 

' In recent years, numerous authors have investigated the phase diagram of symmetrical binary mixtures usin g liquid 
, state theories (of var ying de gree s of sophistication), as well as computer simulations ItLIS 1^ ITol ITl|. IT^ ll^ Il4 
- - ' Hsl ITgL ItR ITsL [1^ I20I I21L 123 . l23 | . However, these studies have been confined for the large part to the special case of 
equimolar concentrations of the two species (i.e., for Afi — 0). The resulting picture of the phase behaviour shows 
four distinct topologies (or archetypes) of phase diagram as 6 is varied. Starting at 5 = 1, and in order of decreasing 
6, these are distinguishable as follows: in the first type, the critical line that characterizes the demixing transition 
(the so-called A-line) intersects the LV coexistence line at liquid densities and at temperatures well below the LV 
critical temperature, forming a critical end point (CEP); in the second type, the A-line approaches LV coexistence 
again from the liquid side, but slightly below the LV critical temperature, terminating in a separate tricritical point 
which itself marks the end of a first order transition between mixed and dcmixcd liquid phases; in the third type the 



2 



A-line intersects the LV critical point directly, the two combining to form a tricritical point. Finally in the fourth type 
the A-line intersects the low density branch of the LV coexistence curve, forming a CEP. For clarity we summarize 
m Table d the various notations that have been introduced in literature to classify the phase diagrams of bi nary 
symmetrical mixtures for A/x = 0. A schematic presentation of the first three types can be found in Fig. 1 of 13|, 
while the fourth type is shown in Fig. 1 of 22] . We note that in analytical calculations the precise value of S marking 
the boundaries between two successive types for a particular interaction potential, depends to some extent on the 
method used for the calculation. 

For the more general case of non-equimolar concentrations (i.e., allowing for A/i ^ 0), relatively few investigations of 
the phase behaviour can be found in literature 0, 0, |3 ■ This fact undoubtedly reflects the obstacles encountered 
by both simulation and theoretical approaches in dealing with the increased complexity of the full phase diagram 
spanned by temperature T, pressure p, and the chemical potential difference A/i. Probably the first study to provide 
a rough idea of the complexity of the full phase diagram was that of Antonyevich et al. > based on conventional 
integral-equation techniques. It was followed shortly afterwards by the remarkable study of Pini et al. who used 
the highly accurate Hierarchical Reference Theory (HRT) [2^, complemented by MF calculations, to provide the first 
insight into the complex topology of the full phase diagram of binary symmetrical mixtures for A/i ^ 0. However, 
the high computational costs of HRT did not permit a systematic discussion of the topology of the phase diagram. 

In this work we present a systematic investigation of the full phase diagram for one of the archetypes of phase 
diagram outlined above, namely that exhibiting both a LV critical point and a first order liquid-liquid transition 
terminating at a tricritical point (27j . Our investigations are based on the mean spherical approximation (MSA) and 
grand canonical Monte Carlo (GCMC) simulations. We have chosen MSA for the simple reason that for a mixture 
of HCY fluids it provides closed (i.e., semi-analytic) expressions for the structural and the thermodynamic properties 
|28j |: these concern in particular the pressure, p, and the chemical potentials, /ii and fi2. The formalism for a general 
mixture of HCYs is summarized in (2^. Using then MATHEMATICA"'"^ [s^I, the solution of the coexistence equations is 
considerably less involved than it is the case for integral-equation theories that have to be solved fully numerically jj] . 
Although MSA is less accurate (in particular in the critical regions) than, e.g., advanced liquid state theories, such 
as HRT 25] or the self-consistent Ornstein-Zernike approximation (SCOZA) 0,0,1^, it is nevertheless able to 
provide results of reasonable accuracy; this will also be corroborated in the present contribution. GCMC, when used 
in combination with extended sampling and histogram reweighting techniques provides a powerful means of studying 
fluid phase equilibria with a precision which is limited only by statistical errors 0, [s^ • We extend its use here to 
the case of a chemical potential difference fleld, thus providing a benchmark against which to test the MSA results. 

Our work makes progress on two fronts. First, we demonstrate that liquid state theories and computer simulations 
have reached levels of accuracy and efficiency that the full phase diagram [i.e., in (T, p, c)- or in (T,p, A/x)-space] of 
a symmetrical binary fluid mixture can be determined on a semi-quantitative level. We note that even a few years 
ago such an aim was considered to be too ambitious |l9l ]. Second, we show that our simple model exhibits a wider 
spectrum of topologies of phase behaviour than one might have expected from previous investigations, in particular 
from those that were restricted to the equimolar case. Specifically, we confirm that for the archetype of phase 
behaviour considered (whose associated domain of 6 is defined on the basis of the A/i = phase diagram topology), 
two distinct subtypes emerge when one considers A/i 0, as already suggested by Pini et al. |19| on the basis of 
selected isothermal cuts through the phase diagram. Our exposition of the full phase diagram sheds further light on 
the nature of the subtypes, by showing that their topologies are defined by the location of four coexistence surfaces 
(each of which carries a critical line), which intersect in four triple lines. Some of these triple lines do not originate at 
the A-line at equal concentrations. We consider the locus of the triple lines to be the key quantities that distinguish 
the two subtypes of phase behaviour. We discuss in detail the transition scenario between the two subtypes, which 
has not yet been documented in literature and which reveals new and interesting details about the critical behaviour 
of the system. The existence of the two subtypes is confirmed by the results of the GCMC simulation studies which, 
while being less comprehensive than the theoretical ones, demonstrate that simulation can be predictive regarding 
quite subtle features of the phase diagram. 

Although a binary symmetrical mixture might at first sight seem to be a system of purely academic interest, argu- 
ments have been put forward in literature that they are able to reflect features of realistic systems: van Konynenburg 
and Scott Q note that 'the only real binary systems in which such symmetry occurs are mixtures of d, Z-optical iso- 
mers 's^', while Woywod and Schoen 24] flnd that the phase behaviour of ^He-^He mixtures contains many features 
of the topology of the phase diagram of binary symmetrical mixtures. Furthermore, as we shall discuss later, the 
symmetrical mixture has direct relevance to some classes of spin fluid models. 

The paper is organized as follows. In Sec. Ill Al we introduce our model, a symmetrical binary fluid mixture in 
which particles interact via the HCY potential. The MSA integral-equation strategy is then outlined, as is the grand 
canonical simulation methodology. Sec. IIIII represents in turn our MSA and GCMC results for the phase behaviour 
of the system for non-zero chemical potential difference fleld. Finally in Sec. IIVI we attempt to set our results within 
their wider context, and assess the outlook for future work. 
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II. MODEL AND METHODS 



A. The model 

We consider a binary symmetrical mixture where the particles interact via the HCY potential 

^ \ -etjcrl exp[-z(r - a)] r Xr ' ^ ^ 

Here the e^j are the contact values of the potentials, a is the hard-core diameter, and z is the inverse screening length, 
assigned the value z<t — 1.8 in this work. Further the system is characterized by a temperature T [with /? = (/cbT)"^] 
and a number density p = N/V, N being the total number of particles in a volume V . Defining c as the concentration 
of species 1, we introduce partial number densities pi = cp and p2 — (1 — c)p; in the expressions below we use Xi = c 
and X2 = 1 — c. Standard reduced units will be used throughout the paper. 



B. Mean spherical approximation 

The solution of the Ornstcin-Zernike equations 3] along with the MSA closure relation, i.e.. 



c^J{r) = -/3$jj(r) r > a , . 

g,,{r) ^ r<a 

can be carried out for this particular system to a large extent analytically . Here the Cij (r) and the gij (r) are the 
direct correlation functions and the pair distribution functions, respectively. Our solution strategy is based on the 
general formalism presented by Arrieta et al. ^29] for the MSA solution of multi-component HCY systems. These 
authors have derived a set of coupled, non-linear equations for the initially unknown coefhcients that characterize the 
correlation functions Cy(r). These equations require a fully numerical solution (see below); for details we refer the 
interested reader to ref. [23 ■ Once these coefficients are obtained, essentially all thermodynamic properties can be 
calculated via simple algebraic expressions. Below we present only those that we require for the present contribution. 

Within the MSA the thermodynamic properties derived via the energy route are usually more accurate than via 
other routes [s^l- In the following, we list expressions for the excess (over the hard-sphere (HS) reference system) 
pressure and chemical potentials. The corresponding properties of the reference system are given by the Carnahan- 
Mansoori-Starling-Leland (CMSL) equation of state With pE = _pg^^^ + ApE and p^ = p™^^ + A/i, (i = 1, 2), 
these quantities are given by 



^ = ^pJ2x.x,amg.,{<j)r [5° (a)]n + J (4) 



-I- 



1 — 77 / 1—7] 

•2 
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(1 — 77)^ 77(1 — 77)^ 

3a2eiln(l-77) (a3e|)(2 ~ 57/ + 77^) 



7^2 77^(1 — 77)^^ 

2a3e3ln(l-77) 



(5) 
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/3A^l, = - — ^ K., G., - 2 E z^^- [c^J- (0) - ^° (0)] (6) 

3 3 

where — 7r/6(p(T-'), j = 1,2,3; in particular ^3 ~ rj, the packing fraction. Further parameters that appear in 
these relations (such as J, Gy , and Kij) are a priori unknown quantities, which are obtained from the numerical 
solution of the non-linear equations mentioned above. The gijicr) and the g^j{<y) are the values of the pair distribution 
functions of the HCY and of the HS reference system at contact. Finally, Cij{Q) and c°^(0) are the zero-components 
of the Fourier-transform of the direct correlation functions for the HCY and the HS system. Details of the numerical 
solution of the MSA which leads to the full phase diagram will be presented in the Appendix. 



C. Computer simulations 

Many features of the simulation techniques employed in the present study have previously been detailed elsewhere 
[llL 13^ . Accordingly, we confine the description of our methodology to its essentials, except where necessary to 
detail a new aspect. 

GCMC simulations were performed using a standard Metropolis algorithm [l3l l36j . The simulation scheme com- 
prises two types of operations: 

1. Particle insertions and deletions. 

2. Particle identity transformations: 1 ^ 2, 2 ^ 1. 

Since particle positions are sampled implicitly via the random particle transfer step, no additional particle translation 
moves are required. Indeed it is the sampling of density fluctuations which represent the bottleneck for this problem, 
and the incorporation of particle translations does nothing to improve sampling efficiency. 

The simulations were performed using the potential of eq. 1^, truncated at rc = 3.0(7. A mean- field long range 
correction for the effects of this truncation was applied in the usual manner Q . To minimize the overheads associated 
with identification of neighbor interactions, we employed a cell decomposition scheme [s^. This involves partitioning 
the periodic simulation space of volume into cubic cells, each of linear dimension of the interaction range, i.e., 
L/l = 3.0. We chose to study system sizes corresponding to Z = 4, / = 7, and I = 9. The smallest system size was 
used to trace phase boundaries away from critical points; additionally, the larger ones were employed near criticality 
to facilitate crude finite-size scaling estimates of critical point locations. At each state point, equilibration periods of 
up to 2 X 10'^ MCS (Monte Carlo sweeps, where one sweep comprises P attempted transfers and identity swaps) were 
used. Sampling frequencies were 500 MCS for the I = 4 system size rising to 2000 MCS for the I — 9 system. The 
total number of samples drawn for each state point studied was about 10®. 

In this work we explore the parameter space spanned by the variables (/ii, A/i, T, S). To accomplish this, without 
having to perform a very large number of simulations, we employed the histogram reweighting technique 37]. Use 
of this technique permits histograms obtained at one set of model parameters to be reweighted to yield estimates 
corresponding to another set of model parameters. To enable simultaneous reweighting in all four fields {ni,Afi, T, S), 
one must sample the conjugate observables {pi,pi — p2,u,Ud), with pi = Ni/V, p2 = N2/V, u the configurational 
energy density, and Ud that part of u associated with interactions between dissimilar particle species. These quantities 
were accumulated in the form of a list during the simulations, and histograms of their forms were created and 
reweighted via post-processing of the entries in the list. 

Within a region of first order coexistence, a standard grand canonical ensemble scheme is severely hampered by 
the large free energy barrier separating the two coexisting phases. This barrier leads to metastability effects and 
prohibitively long correlation times. To circumvent this difficulty, we have employed the multicanonical preweighting 
method [l^ which enhances sampling of the interfacial configurations of intrinsically low probability. This is achieved 
by incorporating a suitably chosen weight function in the GCMC update probabilities. The weights are subsequently 
'folded out' from the sampled distributions of observables to yield the correct Boltzmann distributed quantities. Use 
of this method permits the direct measurement of the distribution of observables at first order phase transitions, even 
when these distributions span many decades of probability. Details concerning the implementation of the techniques 
can be found in references |,38„ .39J . 

Phase boundaries were traced in the space of {pi,Ap,T) by applying the equal peak weight criterion ^3 to the 
distribution of the appropriate fluctuating order parameter, whether this be the concentration or the overall density. 
Critical point parameters were estimated using a crude version of the finite-size scaling (FSS) analysis described in ref. 
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[33 ■ The analysis involves scanning the parameter space until the observed probability distribution of the fluctuating 
order parameter matched the independently known universal fixed point form appropriate to the Ising universality 
class in the FSS limit jij. 



III. RESULTS 
A. MSA results 

In the following we present the MSA results for the two subtypes of phase diagram, denoted a and that we 
have distinguished. For each subtype we discuss its form first in the space of the physical fields, (T,p, A/z), and 
then in a mixed field-density space, (T, p, c). Although representations in field space are somewhat less conventional 
than those in a mixed space of fields and densities, they do afford a compact display of results because two-phase 
coexistence occurs on surfaces. Taken together both representations provide a full picture of the characteristic features 
of the respective topologies. We note that as a consequence of the symmetry of the underlying model the field space 
representation is symmetric with respect to the plane A/i = 0, while the mixed space representation is symmetric 
with respect to c = 1/2. 



1. The type a phase diagram 

Here we consider the phase behaviour of subtype a. Motivated by the results of Pini et al. j^, we have obtained 
the phase behaviour using MSA for 5 = 0.67. We present first, in Fig. ^ the phase diagram in (T, p, A/i)-space. 
Four coexistence surfaces are distinguishable, 7^", i = 1, ... ,4, along each of which two phases coexist. TZ^ lies in 
the (A/i — 0)-plane and represents the demixing surface. T^-j and TZ'^ are symmetrically related and describe neither 
pure demixing transitions nor pure liquid- vapour (LV) transitions. Finally, along T?," there is (predominantly) LV 
coexistence. The lines of intersection between 7?.2 and with TZ^ are symmetrically related triple lines, and i'^. 
On the high temperature side they terminate in critical end points (CEPs), which are connected by a critical line, 
cf, marking the upper edge of surface TZ'^. This critical line passes through the LV critical point of the field- free 
case, located at A/j, = 0. The surfaces 7?,f , TZ^, and TZ^ thus form a pocket delimited by the critical line cf and the 
triple lines t2 and tg. In this pocket a homogeneous, mixed fiuid is stable as was observed previously in ref. |19| . 
Two further triple lines, t" and t" , can be found in the equimolar plane of TZ^: they are the intersection lines of the 
two symmetrically related surfaces, 7?.2 and TZ'^, with the demixing surface TZf. The fact that these three surfaces 
intersect in a line rather than in a point is a consequence of the symmetry of the system. All triple lines meet in the 
only quadruple point of the system. The equimolar triple line t" intersects at high temperatures the tricritical point 
of the system. At this point three critical lines, c", C2 , and C3 , originate: in the equimolar plane of TZ^ lies the critical 
line c" that characterizes the demixing transition (the so-called A-line) and the two symmetrically related critical lines 
C2 and C3 that specify the mixing-demixing transition, and which tend (for |A/i| — > 00) towards the critical points of 
the respective pure components. If we focus on the isotherms that are marked in Fig. ^ we recover at a quantitative 
level the scenario depicted in Figs. 3(b) to (d) of ref. [l9j . 

Turning next to the mixed field-density representation for subtype a [i.e., in (T, p, c)-space], this is displayed in 
Fig. [3 Four distinct coexistence surfaces, 5f , are observable, corresponding to the respective T^f , i = 1 . . . , 4. In 
this representation, however, the surfaces separate 'two-phase' regions (below) from regions of homogeneous phases 
(above). Again, surfaces intersect along triple lines now represented by three lines (branches). Each surface exhibits a 
critical line, cf . 5f is symmetrical with respect to the (c = l/2)-plane and its critical line, c", is the A-line. On 5f , two 
symmetrical high-density phases (p > 0.5), specified by (p, c) and (p, 1 — c), are in coexistence. The two (symmetrically 
related) coexistence surfaces, 5^ and ^3 , are encountered for low and intermediate densities (p < 0.6); the order 
parameter that characterizes this transition is a linear combination of the density and concentration differences of 
the two coexisting phases. Coming from high temperatures the A-line bifurcates at the tricritical point into two 
critical fines, C2 and (belonging to 5^ and 5^), which pass through minima (i.e., 'double critical points') and 
head for the respective LV critical points of the pure phases. The fourth coexistence surface, S^, is predominantly 
LV-like in character and is located in between S2 and Si^ . Its critical line passes through the LV critical point of the 
field-free mixture. The surface as a whole is delimited by the four intermediate-density branches of the triple lines 
t2 and t3 which form a lens-shaped 'loop'. Tie-lines starting at these enclosing triple lines connect a vapour and a 
liquid phase of approximately equal concentrations with a liquid phase of higher density and different concentration 
[see (p, c)-projection of Fig. [21 where for each triple line a set of tie lines was selected]. The latter states form the 
high-density branch of the triple lines tf and tg, that are located in the 'valley (s)' formed by S2 (or S^) and Sf . In 
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(T, p, c)-space the quadruple point can be localized as follows: its intermediate density representations are located on 
the end points of the lens-shaped 'loop', while the other two, being symmetrically related, are the bifurcation points 
of the triple lines and tf. As the critical line c" of 5" approaches the bounding triple lines and tg , the coexisting 
vapour and liquid phases become critical. Since these are simultaneously in equilibrium with a non-critical phase (the 
so-called spectator phase) located at the end point of the high-density branch of the triple lines, this point is a CEP. 
The region of a homogeneous, mixed fluid at intermediate densities found in |19| | and discussed above in field space, 
is bounded by the three coexistence surfaces 5^ , iSg , and S2 ■ 

We close with a more qualitative discussion of the type a phase diagram which, nevertheless, might be rather 
instructive. In Fig. |31we show the projection of the phase diagram in (T, p, c)-space onto the (T, p)-plane. In this 
representation we have introduced six temperature ranges (labeled A' to F'); they are separated by five isotherms 
which are characterized by particular temperatures of the system. In Fig. ^we show representative isothermal cuts 
taken from each of these temperature ranges in the (p, c)-plane. These figures complement those of Fig. 2 in ref. 
in the sense that they present a more comprehensive sequence of possible scenarios of coexistence regions encountered 
in the type a phase diagram. However, we note that within subtype a, other scenarios in the isothermal cuts can also 
be observed. This is, for instance, the case when the LV critical point lies above the temperature that characterizes 
the minimum in the symmetrically related critical lines C2 and C3 that pass through the equimolar tricritical point. 
Hence, scenarios different to the ones presented here can be encountered. In range A', i.e., at temperatures below the 
quadruple point, regions of homogeneous phases are only found at low densities (extending over all concentrations), 
and in small areas at high densities (with concentrations close to or 1). As indicated by the tie lines, two-phase 
coexistence is generally observed; only at c = 1/2 and low densities and at the symmetrically related corners of the 
homogeneous high density phases does three phase coexistence occur. At the temperature of the quadruple point, an 
infinitesimal amount of a fourth homogeneous phase emerges at c = 1/2 and p ~ 0.45. As we proceed to range B', 
this point extends to a small region where a homogeneous phase exists (see the discussion above); as can be seen from 
the tie lines, its three corners are points of three phase coexistence. Furthermore the other regions of homogeneous 
phases encountered already in range A' grow both in density and concentration. 

By further increasing the temperature we pass through the equimolar LV critical point and reach range C. Now 
a narrow passage has opened up between the homogeneous mixed fluid (c 1/2 and p ^ 0.5) and the homogeneous 
low density phases, that becomes broader with increasing temperature. It is characterized by two critical points 
located on the symmetrically related coexistence regions. At some temperature these become CEPs and we enter the 
range D'. Here the scenario is dominated by three huge, interconnected coexistence regions: one is symmetrical with 
respect to c = 1/2 at high density and two are symmetrically related (one for c > 0.6, the other one for c < 0.4). 
In the range D', the coexistence regions are nearly entirely two-phase in character; three-phase coexistence is only 
observed at intermediate and high densities. The next transition temperature is characterized by the minimum of 
the symmetrically related critical lines C2 and C3 that originate at the tricritical point and head towards the critical 
points of the pure components. As we pass this temperature, we enter range E' and the connections between the three 
coexistence regions break, leaving them entirely disconnected. Four critical points can now be observed: two of them 
are located on the symmetrically related coexistence regions and two others (symmetrically related), can be found 
on the demixing coexistence region. In small areas, close to these regions, three-phase coexistence is observed. As 
we finally pass through the tricritical point of the system we enter range F', where the three coexistence regions are 
well-separated and the phase diagram is dominated by huge areas of homogeneous mixtures. The coexistence regions 
now show exclusively two-phase coexistence. The two critical points of the demixing coexistence region have merged 
at the tricritical point to a single critical point (located on the A-line). 



2. The type f3 phase diagram 

In the type (3 phase diagram, which we have found to occur at (5 = 0.69, the phase behaviour is distinct from that 
of type a. In the {T,p, A/i)-space representation of Fig. Owe can again clearly identify four coexistence surfaces 7?.f , 
i — 1, . . . , 4. However, the manner of their intersection engenders a different topology. The triple line is now the 
intersection of the surfaces TZi and TZ^. Intersection of 7^.2 and TZ^ with TZ^ leads again to the triple line t^, in the 
(A^ — 0)-plane, which terminates at high temperatures at the tricritical point. From there, (and similarly to type 
a), three critical lines originate: the A-line (i.e., the critical line cf of TZi ) continues to higher temperatures, while 
two symmetrically related critical lines, and Cg, head off to lower temperatures. These terminate in two CEPs 
where they also meet two symmetrically related triple lines, and tg. Along these lines 7^2 TZ^ intersect with 
7Z^ (which is a predominantly LV coexistence surface) and which dominates the low pressure region and covers the 
whole A/i-range. The intersection of t2 and tg with the aforementioned triple lines and t^, in the (Ap = 0)-plane, 
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gives rise to the only quadruple point of the system. TZ'l is delimited at high temperatures by a further critical line, 
C4, which connects the critical points of the pure components (i.e., for |A/i| — > 00) with the LV critical point of the 
mixture at Afj. — 0. 

Again considering the mixed field-density representation in (T, p, c)-space, the distinction between type P and type 
a phase diagrams is clear. The type (3 case is depicted in Fig. Four coexistence surfaces 5f with accompanying 
critical lines cf can be identified; they correspond to the respective 7?,f . The A-line (cf ), being the critical line of the 
symmetrical demixing surface 5f , bifurcates at the tricritical point into two critical lines, C2 and Cg, which traverse 
the (symmetrically related) surfaces S2 and . Passing through minima (i.e., 'double critical points') these critical 
lines terminate in CEPs located at the symmetrically related high density branches of the triple lines, and t^. These 
triple lines are the intersection of S2 and with coexistence surface iS^* that encompasses the entire concentration 
range. At these CEPs, two high-density phases become critical while in coexistence with the spectator phase located 
at the end points of the low-density branches of the triple lines. The critical line cf on iSf passes through the equimolar 
LV critical point and connects the critical points of the respective pure components; in contrast to type a it is now 
completely detached from the A-line. In subtype f3 the triple lines also show a distinctively different behaviour than 
their equivalents in subtype a. As can be seen from the representative tie lines (all connecting coexistence points on 
triple lines) shown in Fig. El a single low-density phase (with c ~ 1/2) is connected via tie lines to two high-density 
phases. In contrast to type a, the high and the intermediate density triple lines, tf , i^, and tg, merge to form two 
symmetrically related 'loops', while the low density branches terminate as the spectator phase of the CEPs. Again, for 
p ~ 0.5 a region of a homogeneous, mixed fluid is encountered, which is enclosed by the coexistence surfaces S2 , , 
and (see also ref. ^3)- ^hc quadruple point, the low density branch of tf bifurcates in the low density branches 
of and tg , whereas the high density branches of bifurcate in the high density branches of (and symmetrically 
related tg) and tf . The fourth representation of the quadruple point is located in the 'valley' of mixed fluid, where 
the intermediate branches of , tg , and tf intersect. 

Again we close the discussion of the type f3 phase diagram with a more qualitative discussion. In Fig. [7|we show 
projections of characteristic lines of the phase diagram onto the (T, p)-plane, defining six temperature ranges, labeled 
A to F. In Fig. |H1 we show isothermal cuts for each of these temperature ranges. In the type f3 phase diagram 
the scenarios in the temperature ranges A and B (which are again separated by the quadruple point temperature) 
correspond on a qualitative level exactly to the corresponding ones in type a. In particular, in range B we again 
encounter three large interconnected coexistence regions with a small area of a homogeneous, mixed fluid for c ^ 1/2 
and p ~ 0.5. The minimum on the symmetrically related critical lines, and Cg, on S2 and iSg (mentioned above) 
define the temperature where we enter range C: in contrast to type a, the three large coexistence regions break up 
at p > 0.5, splitting off the demixing coexistence region and leaving the other two connected. The symmetrically 
related passages are characterized by pairs of critical points (located on each of the two coexistence regions in a 3d 
representation in (T, p, c)-space) , and connected by symmetrically related critical lines C2 and Cg that originate in the 
tricritical point and terminate at symmetrically related CEPs. This CEP-temperature marks the limit between ranges 
C and D. In range D only the symmetrically related critical points on the demixing surface remain; simultaneously, the 
connecting bridge between the large coexistence regions that extend over all concentrations becomes narrower. Passing 
the equimolar LV critical point (and thus entering range E) this connection breaks up leading to two symmetrically 
related coexistence regions that are entirely detached from the demixing region, and on each of which a pair of 
symmetrically related critical points emerges. Finally, above the tricritical point of the system (i.e., in the range F) 
the three coexistence regions are by now well separated, each being characterized by a critical point, which corresponds 
exactly to the scenario depicted for the range F' in subtype a. Again we note that within subtype slightly different 
scenarios in the isothermal cuts can be observed. This is, for instance, the case if the equimolar LV critical temperature 
is lower than the temperature of the minimum in the symmetrically related triple lines, and tg. 

3. The transition scenario 

Of course it is of particular interest how the transition from one subtype to the other takes place. We start from 
subtype a and increase the parameter 6. In the (T,p, A/^)-representation the transition scenario between the two 
subtypes can be traced easily. The characteristic pocket formed by 7^2 j ' ^"^^ ^4 becomes larger by extending the 
limiting triple line branches, t2 and tg (and thus the connecting critical line C4 ) in both directions. Simultaneously the 
symmetrically related critical lines C2 and Cg that emerge from the tricritical point are shifted to lower temperatures. 
At the transition value S, these critical lines have reached the CEPs, thereby forming two additional, symmetrically 
related tricritical points. As a consequence the coexistence surface 72.2 is subdivided by the triple line \^ (the former 
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tj), into two sub-surfaces, 7^2^ and Ti}^"- The same happens with the symmetrically related surface 7^" . A three sided 
pocket is now formed by T^-j//, T^y^^ and TV^ (the former TZ1). From below it is bounded by the quadruple point, its 
edges are the triple lines, the upper three corners are the tricritical points which, in turn, are connected by sections 
of the critical lines (see below). Upon further increasing 5, TZf merges with 7^2^ and TZ^^, forming thus TZ^ while 7?.2^/ 
(TZyi) becomes T^f (^a)- The surface TZf is not affected by the transition and corresponds to TZi . In a similar 
manner the critical line C2 (and symmetrically related C3 ) is subdivided at the transition in two sections d§ and C2rr , 
where is delimited by the field- free tricritical point and one of the newly formed tricritical points, c^"^ is limited by 
the newly formed tricritical points. Upon further increasing i5, C4'', c^,,, and Cy, merge to form cf. (cy) becomes 
C2 (C3) and are delimited by CEPs and the field-free tricritical point. 

To discuss the transition scenario in the (T, p, c)-representation we again start from type a and increase 6. Then 
gradually the high- and the intermediate-density branches of the triple lines {2 and t" lengthen until they merge when 
the CEP located on the intermediate-density branch meets the end point of the high-density branch, thereby forming 
a tricritical point (see Fig. |5J). Upon further increase of S, a high-density branch (with a CEP) and a low-density 
branch of the triple line detach, thus resulting in the topology of the /3-subtype. Concomitant with this metamorphosis 
of the triple lines, we observe a related development in the critical line c" that passes through the LV critical point 
in subtype a: with increasing 5, this line lengthens in both directions until it meets (at the crossover between the 
subtypes) the critical lines of the surfaces S2 and S^. As S is further increased, a new critical line forms connecting 
the critical points of the pure phases, passing through the equimolar LV critical point so that the CEPs remain with 
the newly formed triple lines at higher densities. 

Within the MSA the transition occurs for Smsa = 0.678(0). We derive this estimate by following the development 
of the triple lines t2 and to {2 and t'^ with increasing 5. In both subtypes, two of the three branches of the triple 
lines mentioned above meet each other in CEPs, whereas at the transition all three branches meet in a tricritical 
point. Fig. shows a selection of the triple lines t3 and tg for different values of S. 

The transition scenario (although not explicitly identified there) is nicely depicted on several occasions in [l^ . 
Within the MF framework, where the transition occurs at Smf = 0.65338, the above mentioned loop of critical lines is 
nicely depicted in Fig. 1(c). Figs. 10 and 12 show projections of the HRT (T, p, c) phase diagram onto the {p, c)-plane 
at various temperatures for 6 — 0.67, while in Fig. 11 isothermal cuts through the (T,p,Ap) phase diagram are 
shown. This S-vahie is close to the HRT-transition value, (5hrt- 



4. Relationship to previous theoretical work 

The two subtypes that we discuss here have already been classified by van Konynenburg and Scott type a 
corresponds to 'sym. Ill- A*' while type (3 to type 'sym. II-A*'. Several years later Antonevych et al. [lj| pointed 
out that the complex topology of the full phase diagram of a binary symmetrical mixture originates from the complex 
interplay of the LV transition and the demixing transition (see Figs. 1 to 4 in this contribution). However, most of 
their considerations were carried out on a qualitative level (see section III.B of ref. |0|). 

Shortly afterwards, Pini et al. presented their remarkable HRT and MF theoretical study of the full phase 
diagram of binary symmetrical mixtures, discovering several interesting features which they presented in the form 
of isothermal cuts through the phase diagrams in (T, p, A/i)- and in (T, p, c)-space. The two subtypes that we have 
discussed here in detailed were identified in their contribution as the two alternative ways in which the region of the 
homogeneous fluid occurring at intermediate densities (see above) is formed from the intersection of the coexistence 
surfaces. Additionally, the authors pointed out that in what we term subtype a, the critical temperature at A/i — 
is always lower than the tricritical temperature, while in subtype f3 no such ordering relation holds. While Pini et 
al. focus in their contribution on a classification of the topology via the critical lines, we consider a discussion of the 
phase diagrams via the location of the triple lines to be more enlightening. In fact, our 3d representations both in 
the space of the physical fields as well as in the mixed field-density space helped us to clarify the following points: 
(i) from these representations we obtained an unambiguous identification of two distinct subtypes (confirming also 
the conclusions drawn in J,9]); (ii) in the 3d plots, the 'double critical points' reported in J^] are identifiable as local 
minima in the critical lines; (iii) we find previously unreported triple lines; (iv) the transition scenario between the 
two subtypes could be discussed and clarified in detail (see sec. IIII A3|) . Comparison of our isothermal cuts with 
those presented by Pini et al. facilitates the link between our subtypes and the topological changes identified in |l9|: 
in their Figs. 2, 3, 6, and 7 a type a behaviour can be identified, while Figs. 8 and 9 display a (3 topology; Figs. 10 
to 12 depict a close-to-transition scenario (see also sec. IIII A^ . 

Finally, Woywod and Schoen 24] have very recently studied a lattice model using a MF density functional theory, 
to calculate the complete phase diagram of a binary fluid mixture of equally sized particles. The conclusions that the 
authors could draw from their investigations are relevant also for the present study: they present arguments (that back 
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up earlier conclusions |42|) which preclude the existence of tricritical points in general binary mixtures, but allow them 
under certain symmetry conditions which are fulfilled by symmetrical mixtures, but can also be found in particular 
asymmetric mixtures such as those investigated in 12^. Their investigations of the binary symmetrical mixture is 
focused on the shape of the phase diagram in field space for selected values of a parameter A, that characterizes the 
relative strength of the mixed interaction (cf. thus following essentially the sequence of systems summarized in 
Tabled However, the existence of the two subtypes a and /? is not mentioned at all in this contribution. However, due 
to the fact that Woywod and Schoen could calculate thermodynamic potentials via an essentially closed expression, 
they were able to pin down exactly the transition scenario discussed in sec. IHI A 31 In particular in their Fig. 5 
this transition in the (T, /i, A/i)-space is clearly depicted, showing explicitly the three tricritical points, the quadruple 
point, and the connecting critical and triple lines that occur for this particular system. We point out that for the 
semi-analytical approach used in the present contribution within the MSA framework, such an exact localization of 
the transition is beyond reach. 



B. Simulation results 

Turning now to the results of our GCMC simulations, our strategy for mapping the coexistence behaviour was 
initiated on the symmetry plane A/i = 0, where we first mapped the phase diagram in the (/ii, T)-plane. Use of 
histogram reweighting then provided estimates for the phase behaviour at small, but finite, Afi. Guided by this 
prediction, a new set of simulations were subsequently performed at near coexistence state points for this value of A/i, 
the results of which were extrapolated to yet larger A/i, and so on. In this manner we were able to track the phase 
behaviour as a function of Afi and cover a large range of concentrations in the process. By accumulating separately 
contributions to the energy from like and unlike particle interactions, we were further able to perform histogram 
extrapolation with respect to S. This was useful in helping to find the regions of S relevant to the two subtypes. 

The high dimensionality of the full phase diagram precludes a comprehensive study of the full 3d phase diagram 
as was done via MSA. Nevertheless projections of our data onto the (T, /?)-, (p, c)-, and (T, c)-planes, as shown in 
Figs. El and Figs. ^1 clearly confirm the existence of the two subtypes seen in the corresponding projections of 
the MSA data (Fig. |21 and EJ and detailed above, albeit at slightly different values of S. Specifically, for S = 0.66, 
the critical lines emanating from the tricritical point links up with that coming from the LV critical points of the 
pure phases, while the critical line emanating from the field-free LV critical point terminates at two symmetrically 
related CEPs at which a demixcd high density non-critical (spectator) phase coexists with a lower density critical 
phase. The phase that coexist at the CEPs are joined in Figs. llUf aV^cl via a dashed line. The character of the 
CEP behaviour is clarified by the measured forms of the density and concentration distributions, as shown in Fig. 1111 
The form of P{p) for the critical phase has a highly non-Gaussian structure which matches the universal Ising order 
parameter distribution appropriate for a critical finite-sized system (see inset) [Til Il6^. '32], That the critical phase is 
predominantly liquid- vapor in character is evidenced by the absence of strong concentrations fluctuations as reflected 
in the near-Gaussian form of P{c) for the lower density phase. 

For S = 0.68, the results of Figs. 1121 shows that the phase behaviour is of subtype (3. Speciflcally, the LV critical 
point of the pure phases joins smoothly to that of the equimolar mixture, while the critical line emanating from the 
tricritical point terminate at two symmetrically related CEPs, the coexisting phase of which are joined by dashed lines 
in Figs. El The nature of these CEPs is again elucidated by the corresponding forms of P{p) and P{c) (see Fig. I13|l . 
One flnds now that the spectator phase is a low density mixed phase, while the higher density critical phase exhibits 
fluctuations that are neither predominantly density-like nor concentration-like in character, but a coupled mixture of 
the two. Thus the forms of P{p) and P{c) for the critical phase both exhibit forms that are highly non-Gaussian in 
nature, though neither matches well the universal Ising form. In such situations one can expect that the true order 
parameter for this transition involves a linear combination of density, concentration (and energy), though we have 
not attempted to investigate this matter further here. 

As regards the range of S in which the two subtypes of phase behaviour occur, the simulation results are in semi- 
quantitative agreement with the MSA calculations. However, owing to the computational burden of this problem, 
we were not able to pin down precisely the value of S at which the transition between subtypes occurs, although our 
results shows that 5 differs from that found in MSA by at most 0.01. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have investigated the full phase behaviour of a symmetrical binary fluid in a range of model 
parameters for which, at equal chemical potentials, the system exhibits a LV critical point and a tricritical point. The 
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phase behaviour is considerably richer in both variety and character than would be expected on the basis of knowledge 
of the field-free case alone. Our results confirm previous less detailed reports that for unequal chemical potentials, 
two subtypes a and /3 of phase behaviour occur, and we have elucidated the differences between these subtypes in 
terms of the topologies of triple lines and critical lines. 

The MSA results are in semi-quantitative agreement with those of MC simulations, demonstrating that MSA pro- 
vides a correct description of this system. Use of the MSA should therefore prove useful in narrowing down simulating 
searches of parameter space when seeking a given type of phase behaviour. Moreover, our study has demonstrated 
that simulations are competitive with theory in providing (within a reasonable amount of time) precise information 
on intricate phase diagrams exhibiting complex topologies of critical lines. Given that the various commonly used 
theoretical approaches (HRT, SCOZA, MSA, . . .) do not always agree, MC thus provides an invaluable benchmark 
with which to compare. We note that our results are in good agreement with the less comprehensive MF and HRT 
study of these two subtypes in ref. Although HRT is generally more accurate than MSA, it is very laborious 

to implement and computationally expensive, and still does not produce results in fully quantitative agreement with 
simulation, as has been observed in studies of the field- free case ^^ . 

As regards the more general relevance of our findings, one can consider the symmetrical binary fluid model as 
probably the simplest member of a class of one component 'spin'-fluid models in which the particles carry an orien- 
tational degree of freedom which features in their interparticle potentials |l3l Il9||. Other, more complex examples 
include Heisenberg fS] , Ising , or dipolar spin fluids jl^ . It is well known that similar sequences of phase diagram 
topologies can arise in all members of this class and that one-to-one correspondences can generally be made between 
the phases of the binary symmetrical mixture and those exhibited by the more complex models. However studies of 
phase behaviour for the more complex class members have to be performed using more complicated techniques than 
the MSA. Thus the symmetrical mixture plays a key (computationally tractable) role in elucidating generic aspects 
of the phase behaviour. 

As regards future work, it would be of interest to extend the present studies to encompass the more general case of 
asymmetrical mixtures. Steps in this direction, have recently been reported for the case of a lattice based binary fluid 
model 11^, though this is unable to represent particle species of unequal sizes, whose packing effects are likely to be 
subtle. From the simulation point of view, GCMC simulations of asymmetrical mixtures are not significantly more 
challenging than the symmetrical case, provided the size asymmetry is moderate. Recent methods for minimizing 
finite-size effects in GCMC measurements of coexistence properties of fluid mixtures, should also help render this a 
practical proposition [4^ . 



APPENDIX: NUMERICAL SOLUTION OF THE MSA 

We were able to perform all our MSA calculations using MATHEMATICA"'"^ [30| in a rather straightforward fashion. 
Phase diagrams of general (i.e., non-symmetrical) binary mixture are also accessible via this numerical route. 

We start with the determination of the phase diagram in (T, p, c)-space. To determine at a given temperature T, 
the concentrations and densities of two coexisting phases (labeled I and H) one solves the following (coexistence) 
equations numerically, e.g. with a Newton-Raphson algorithm 

A/if-")(z) = (A.l) 
A/i^")(z) = (A.2) 
Ap(i-")(z) = 0, (A.3) 

here z — {ci, cn, pi, pn; T}. The Afif i = 1, 2, are the differences between the chemical potentials of the respective 
species of the coexisting phases and Ap'^^^^' is the difference between the corresponding pressures. 

At fixed T, equations ljA.l|) to (jA.3|) represent three relations between four unknown quantities, i.e., ci,cii,pi, and 
pii. To solve this set of nonlinear equations one requires an additional constraint, which can contain any subset of 
the unknown quantities. The proper choice of such a constraint is important since it has a distinct influence on the 
efficiency of the calculation of coexistence curves. The simplest (and most obvious) restriction is to fix the value of a 
particular variable and use it to parameterize the coexistence curve. For example, if we wish to calculate, at a given 
T, a LV coexistence line, we may fix the value of the concentration of the vapour phase, ci, and solve equations HA.1() 
to ljA.3|l for cii, Pi, and pu. In order to proceed to a neighboring pair of coexisting states, we then add to this solution, 
z say, a set of parameters Az = {Ac, 0, 0, 0; 0}, i.e., z' = z + Az, which constitutes a starting point for solving the 
coexistence equations for the new pair of coexisting states. The step size Ac in the concentration ci depends on the 
form of the coexistence lines. Solution of the coexistence equations leads to a neighboring pair of coexisting states. 
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Given any two states, new initial values for coexistence points can be obtained via linear extrapolation from the 
previous states, thus facilitating use of a larger step-size in Ac. 

For a triple line (i.e., for three-point coexistence) we have to satisfy, in addition to equations l)A.l|l to (|A.3|I . the 
relations 

A^f^""^'(z) = (A.4) 
AmP-"^'(z) = 
Ap("-"i)(z) ^ 

(A.5) 

with z being now z = {cj, cn, cni, Pi, Pii, Piii; 7"}- For fixed temperature T we thus have six equations for six unknowns, 
i.e., no additional constraint is required. 

The above considerations provide a recipe for locating coexisting states, once a solution is known. An appropriate 
starting point for the iteration outlined above is a coexistence point of the pure fluids, i.e., for c = or c = 1, for 
which either equation (|A.1|I or ljA.2|) is trivially satisfied. Starting at a rather low temperature coexisting state of the 
pure fluid, we employ the above stepwise approach to calculate a full isothermal cut. While for these temperatures 
the coexistence lines of the vapour phase reach the symmetric plane at c = 1/2, this is not the case for the coexistence 
line of the fluid phase: they terminate instead at the high-density branch of the symmetrically related triple lines 
ii and ta. By gradually increasing the temperature we can trace this pair of triple lines to higher temperatures and 
use them at the same time to determine the isothermal cuts of the demixing coexistence surface, S\. For this region 
of the phase diagram, we again take advantage of the symmetry of the system. The variable z is then given by 
z — {cii,cni(= 1 — Cn), Pii, Piii(= Pii)}- In this case equations ljA.l|) and \K.1\ are identical and relation ljA.3|) is 
trivially fulfilled for equal densities. Using coexisting states at a given density and temperature as starting points, we 
scan the demixing surface by gradually increasing the temperature at fixed density. 

Particular care has to be paid to the determination of the critical points which, as a function of temperature, form 
the critical lines. Since MSA can become numerically unstable in the vicinity of a critical point, we have sought to 
localize critical lines within narrow intervals. To this end we have calculated coexistence states characterized by a 
fixed, narrow 'distance' d, with = (ci — Cn)^ -|- (pi — pn)^; for d we have assumed a value of 0.04. In this manner a 
good estimate for the critical line is obtained from the mean value of the concentrations and densities of the coexisting 
phases that fulfill the constraint. 

The procedure outlined above permits construction of the four coexistence surfaces, SJ (7 = a or /3 and i = 1, . . . , 4). 
The final shape of the phase diagram is then obtained from the intersections of these coexistence surfaces; this leads 
- as we have learned from the discussion above ~ to the triple lines. These triple lines, which are a key feature 
of the phase diagram, constitute the boundaries of the coexistence surfaces, and truncate metastable regions of the 
coexisting surfaces, thus defining the final shape of the full phase diagram. 

Once the phase diagram is determined in (T, p, c)-space, we calculate for coexisting states the corresponding ther- 
modynamic properties (notably the pressure and the chemical potential) and determine with these quantities the 
phase diagram in (T,p, A/i)-space. 
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van Konynenburg and Scott 


Tavares et al. 


A < 


T A 

i-A 


Of > 1 


no demixing 


yv > u 


TT A 


a <i 


T 

1 




II-A* 




II-Q 




III-A* 








III-HA 




III 




not classified 




IV 



TABLE I: Correspondence between the types of binary symmetric mixtures classified by van Konynenburg and Scott (see Figs. 
1 and 38 of 1]) and the types I, II, III, and IV introduced by Tavares et al. 4_G|. Types Il-a and II- /3 are subtypes of type II that 
are described in more detail in j4^. A and a are parameters that characterize the relative strength of the mixed interaction of 
the binary (size-)symmetrical mixture. Note that this table was already published in [2^. 
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FIG. 1: MSA results for the type a phase diagram of the binary symmetrical mixture {5 = 0.67) considered in this study in 

(T, p, A/i)-space. Symbols: thin full lines - isothermal coexistence lines, grey full thick line - critical line cj passing through 
the LV critical point, grey dashed thick lines - critical lines (ci , C2 , and cf) passing through the tricritical point, black thick 
lines - triple lines t", i = 1, . . . , 4. 
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FIG. 2: MSA results for the type a phase diagram of the binary symmetrical mixture (5 = 0.67) considered in the study in 
(T, p, c)-space and its projection onto the (p, c)-plane. Symbols: thin full lines - isothermal coexistence lines, dashed thin lines 

- tie lines, grey full thick line - critical line cj passing through the LV critical point of field-free case, grey dashed thick lines 

- critical lines (cf , C2 , and cf) passing through the tricritical point of the field-free case, black thick lines - triple lines t", 
i = 1, . . . ,4. 




FIG. 3: MSA results for the type a phase diagram of the binary symmetrical mixture [5 = 0.67) considered in the study, 
projected onto the (T, p)-plane. Dashed thin lines separate ranges A' to F' (see text). Line symbols see Figs. Q and El in 
addition; thin line - azeotropic (i.e., field-free coexistence) line. 
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FIG. 4: Six different panels display MSA results for the type a phase diagram of the binary symmetrical mixture (5 = 0.67) 
considered in the study, projected onto the (p, c)-plane. The six panels labeled A' to F' show isothermal cuts at the temperatures 
T = {!., 1.035, 1.041, 1.06, 1.07, 1.15} and correspond to the six temperature ranges indicated in Fig. |3 




FIG. 5: MSA results for the type /3 phase diagram of the binary symmetrical mixture {5 = 0.69) considered in this study in 
(T, p, A/i)-space. Symbols: thin full lines - isothermal coexistence lines, grey full thick line - critical line passing through 
the LV critical point, grey dashed thick lines - critical lines (cf , C2, and c§) passing through the tricritical point, black thick 
lines - triple lines tf , i = 1, . . . , 4. 
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FIG. 6: MSA results for the type (3 phase diagram of the binary symmetrical mixture (5 = 0.69) considered in the study in 
(T, p, c)-space and its projection onto the (p, c)-plane. Symbols: thin full lines - isothermal coexistence lines, dashed thin lines 

- tie lines, grey full thick line - critical line cf passing through the LV critical point of field-free case, grey dashed thick lines 

- critical lines (cf , , and Cj) passing through the tricritical point of the field-free case, black thick lines - triple lines tf , 
i = l,...,4. 
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FIG. 7: MSA results for the type P phase diagram of the binary symmetrical mixture {5 — 0.69) considered in the study 
projected onto the (T, p)-plane. Dashed thin lines separate ranges A to F (see text). Line symbols see Figs. |S]and|S] in 
addition: thin line - azeotropic (i.e., field-free coexistence) line. 




FIG. 8: Six different panels display MSA results for the type (3 phase diagram of the binary symmetrical mixture {S = 0.69) 
considered in the study, projected onto the (p, c)-plane. The six panels labeled A to F show isothermal cuts at the temperatures 
T = {1, 1.03, 1.036, 1.05, 1.07, 1.15} and correspond to the six temperature ranges indicated in Fig. |7| 
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FIG. 9: View of the triple lines tf and tf from the equimolar plane (c = 1/2) in (T, p, c)-space for seven different values of 
S = {0.67, 0.675, 0.677, 0.678, 0.6785, 0.68, 0.69} denoted in gray scales from light gray to black. 




FIG. 10: Projections of the phase diagram measured by simulation for S = 0.66 (a) (T, c)-planc. (b) (T, p)-planc (c) (c, /5)-plane. 
Solid lines are first order phase boundaries. Thin solid coexistence lines are influenced by flnitc-sizc effects and serve as a guide 
to the eye. Diamonds represent points on critical lines passing through the field-free tricritical point, whereas circles represent 
points of the critical line passing through the field- free LV critical point. Critical points that belong to the same critical line are 
connected via dotted lines. Large filled symbols represent CEPs and small filled symbols the corresponding spectator phases. 
Tie lines that connect the spectator phase with the CEP are shown as dashed lines. 
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FIG. 11: Critical end point data for S = 0.66. (a) the density distribution P{p); also shown (inset) is the accord between P{p) 
for the critical phase and the appropriately scaled universal Ising fixed point form (see text) . (b) the concentration distribution 
P(c). 




FIG. 12: Projections of the measured phase diagram for 5 — 0.68 (a) (T, c)-plane, (b) (T, p)-plane, and (c) (c, p)-plane. For an 
explanation of symbols, see Fig. 1101 
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FIG. 13: Critical end point data for S = 0.68. (a) the density distribution P(p); also shown (inset) is the accord between P{p) 
for the critical phase and the appropriately scaled universal Ising fixed point form (see text), (b) the concentration distribution 
P(c), where for clarity the two peaks are presented on different scales. 



